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Abstract 



In this paper we prove that near the equihbirum position any periodic 
FPU chain with an odd number of particles admits a BirkhofF normal 
form up to order 4, and we obtain an exphcit formula of the Hessian of 
its Hamiltonian at the fixed point0 

1 Introduction 

In this paper we consider FPU chains with N particles of equal mass (normalized 
to be one). Such chains have been introduced by Fermi, Pasta, and Ulam [5], 
as models to test numerically the phenomenon of thermalization as the number 
of particles gets larger and larger. A FPU chain consists of a string of particles 
moving on the line or the circle interacting only with their nearest neighbors 
through nonlinear springs. Its Hamiltonian is given by 



where y : M ^ M is a smooth potential. The corresponding Hamiltonian 
equations read (1 < n. < A^) 



Here g„ denotes the position of the n'th particle relative to its equilibrium 
position, Pn is its momentum, and throughout this paper we assume periodic 
boundary conditions 
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{qi+N,Pi+N) = {qi,Pi) Vi G {0, 1}. 
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1 INTRODUCTION 



Without loss of generality, the potential y : R ^ R is assumed to have a 
Taylor expansion at of the form 

V{x) = . (ix^ + + + . . .) , (2) 

where k is the (linear) spring constant normalized to be 1 and a, /? e R are 
parameters measuring the strength of the nonlinear interaction. Substituting 
the expression ^ for V into ([1]), the corresponding expansion of Hy is given 

by 

n—1 11—1 71—1 71—1 

(3) 

For any FPU chain, the total momentum P — X]n=iPn integral of 

motion, 

N N 
n—1 n—1 

and therefore the center of mass Q — jj '^n=i In evolves with constant velocity. 
Hence any FPU chain can be reduced to a family of Hamiltonian systems of 
2N — 2 degrees of freedom, parametrized by the vector of initial conditions 
(Q, P) E R^ with Hamiltonian independent of Q. In particular, for N = 2 any 
FPU chain is integrable. Further note that for any vector (Q, P) e R^, the origin 
in R^-2 is an equilibrium point of the reduced system. The momentum of such 
an equilibrium point is given by the constant vector (pi, . . . ,pn) — P (1, . . . ,1). 

Introduce the functions / = (/fe)i<fc<Ar-i, J = {Jk)i<k<N-i, and M — 
{Mk)i<k<N-i defined on R^^^^ with values in R^^-'^ given by 

Ik = ^(xl+yl); Jk = ^(xkXN-k+ykyN-k); Mk = ]^{xkyN-k-XN^yk)- (4) 
Further define the function Ha : R^-i R, given by 

; 1 ^-1 a 2 7 

/r\ ^ . KIT ^ 1 ^-^ p ~ a" s-^ . iTT . TOTT ^ ^ , , 

HaAI) ■■= 2 + 47V ^ + -JlT ^ '''' N aT^'^™' 

k=l k=l '5^™ 

l<i,m<JV-l 

where Ck = Ck{a, (3) := + {/3 — a^) sin^ 
The main result of this paper is 

Theorem 1.1. If N > 3 is odd, any periodic FPU chain admits a Birkhoff 
normal form of order 4 (included). More precisely, for any odd N > 3 there 
are canonical coordinates {xk,yk)i<k<N-i so that the Hamiltonian of any FPU 
chain, when expressed in these coordinates, takes the form 

^+HaAI) + Oi\ix,y)\') 

with Ha.p{I) given by 
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Corollary 1.2. Near the equilibrium state any FPU chain with an odd number 
N of particles can be approximated up to order 4 relative to its center of mass 
coordinates by an integrable system of N ~ I harmonic oscillators which are 
coupled at fourth order. 

Denote by Qa,i3 the Hessian of Ha.pil) at / = 0. Note that Qa,p is an 
[N — 1) X (TV — 1) matrix which only depends on the paramters a and [3. For 
the following result we do not have to assume that N is odd. 

Theorem 1.3. (i) For any given a G K \ {0}, det{Qa^i3) is a polynomial in 
(3 of degree N — 1 and has N — 1 real zeroes (counted with multiplicities). 
When listed in increasing order the zeroes (3k = /3fc(a)(l < k < N — 1) 
satisfy 

0<(3i<a^ <(32<...< I3n-i 
and contain the l^^^^j distinct numbers 

(^1 + (sin^ — )-i j (l</c<L^— j). 

Moreover, index{Qa,i3), defined as the number of negative eigenvalues of 
Qa.,p, is given by 

r 1 for (3 < Pi 

index (Qa.fj) = N for (3i < (3 < (32 

[n^2 for(3> (3n-i 

(a) For a = 0, det((3o./3) is a polynomial in (3 of degree N — 1, and (3 — 
is the only zero o/ det((5o./9)- H has multiplicity N — 1, and the index of 
Qo,f3 is given by 

. , , „ , r 1 /or /3 < 

»nfc(Qo,/3) = |^_2 forp>0 

FPU chains with an even number of particles typically (i.e. ii (3 ^ a^) do 
not admit a Birkhoff normal form up to order 4 due to resonances. Our analysis 
of odd FPU chains leads in the case of even FPU chains to a resonant Birkhoff 
normal form up to order 4 which we use in subsequent work [6J to show that 
their Hamiltonians truncated at fourth order are nevertheless integrable systems 
in the sense of Liouville. 

Recall the definiton ^ of the functions J and M. Let 

Ro^A-lM) := (r{J,M) + Rn{J,M)) (6) 

where 

i?(J, M) = 4 E ^ {jkJM_^ - MkM._,) 
i<fc<f 
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and 

Rn (J, M) 



Jl - Ml if ^ e N 

4 4 

otherwise. 



Theorem 1.4. If N>A is even, there are canonical coordinates {xk,yk)i<k<N-i 
so that the Hamiltonian of any FPU chain, when expressed in these coordinates, 
takes the form 

+ Ho.,p{I) - R^,p{J,M) + 0{\{x,y)\^), 

where Ha^pil) and Ra^p{J, M) are given by @) and (0), respectively. 

Remark 1.5. Theorem \1.4\ can be used to show that a version of Theorem 
\1.1\ holds for N -particle FPU chains (N even or odd) considered with Dirichlet 
boundary conditions by embedding such systems into an invariant submanifold 
of a periodic system with 2N + 2 particles - see JM for details. 



Applications: In the case where N is odd, Theorems 11.11 and 11.31 aUow to 
apply for any given a e M the classical KAM theorem (see e.g. [9]) near the 
equilibrium point to the FPU chain with Hamiltonian Hy for a real analytic 
potential V ^ Ix"^ + + |.t'' + . . . with /3 e R \ {(3i{a), . . .,PN-i{a)}. 
Moreover, as for any given a e M\ {0}, Qa^p is positive definite for (3i{a) < f3 < 
I32{a), one can apply Nekhoroshev's theorem (see e.g. [TOl) to the FPU chain 
with Hamiltonian Hy for V with such /3's. These perturbation results confirm 
long standing conjectures - see e.g. [T|. 

Related work: Theorem 11.11 improves on earlier results of Rink [11] and to- 
gether with Theorem 11.31 solves all open problems stated in [11] for N odd. 
Instead of using symmetry properties of FPU chains employed in |llj our ap- 
proach has been shaped by our earlier work on the Toda lattice [HE]. The latter 
one, introduced by Toda [13] and extensively studied in the sequel, is a special 
FPU chain which is integrable. It turns out that (almost) the same canonical 
transformations which near the equilibrium bring the Toda lattice into Birkhoff 
normal form can be used for any FPU chain. Put in other words, the existence 
of the Birkhoff normal form stated in Theorem ll.il is. at least partially, a con- 
sequence of the fact that the family of FPU chains, parametrized by a, f3, . . ., 
contains an integrable system, namely the Toda lattice. 

Outline: In section |2l we review the notion of a Birkhoff normal form. We 
show Theorem 11.11 in sections |3]|5] and Theorem 11.41 in section |6l whereas Theo- 
rem [US] will be proved in section [71 

Acknowledgement: It is a great pleasure to thank Yves Colin de Verdiere 
and Percy Dcift for valuable comments. 



2 Birkhoff normal form 

Consider an isolated equilibrium of a Hamiltonian system on some 2n-dimensio- 
nal symplectic manifold, i.e. an isolated singular point of the Hamiltonian vec- 
tor field. Neglecting an irrelevant additive constant, the Hamiltonian, when 
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expressed in canonical coordinates w = (g, p) near the equilibrium with coordi- 
nates q = 0, p — 0, then has the form 

= i {Aw, w) + ... 

where A is the symmetric 2n x 2n-Hessian of at and the dots stand for 
terms of higher order in w. We now assume that the equilibrium point w = 
is elliptic, i.e. the spectrum of the linearized system, w = JAw, is purely 
imaginary, spec( JA) — {±iXi, . . . , ±iA„} with real numbers Ai, . . . , A„. Here 

J = ^ ^o" ^ ^^'^ standard symplectic structure of M^". If spec(Jv4) 

is simple there exists a linear symplectic change of coordinates which brings 
the quadratic part of the Hamiltonian into normal form. Denoting the new 
coordinates by the same symbols as the old ones one has 

n 

{AW,W) =J2>'^i1^ +P^)■ 

Definition 2.1. A Hamiltonian H is in Birkhoff normal form up to order 
^ 2, if it is of the form 

H = N2 + Ni + ... + Nm + H,n+i + . . . , (7) 

where the Nk, 2 < k < m, are homogeneous polynomials of order k, which are 
actually functions of qf + Pi, . . . , p^, and where Hm+i + . . . stands for 

arbitrary terms of order strictly greater than m. If this holds for any m, the 
Hamiltonian is simply said to be in Birkhoff normal form. 

Note that if a Hamiltonian H admits a Birkhoff normal form of order m, 
the coefficients of the expansion ([7]) up to order m are uniquely determined, as 
long as the normalizing transformation is of the form id + . . . . However, the 
normalizing transformation is by no means uniquely determined. 

There are well known theorems guaranteeing the existence of a Birkhoff 
normal form up to order m assuming that the frequencies Ai,...,A„ satisfy 
certain nonresonance conditions - see e.g. Theorem 4.3 in [7]. We do not state 
these theorems here, because we will show the existence of a Birkhoff normal 
form up to order 4 (in the case where N is odd) of any periodic FPU chain by 
explicit calculations. 



3 Relative coordinates 

We start by expressing the FPU Hamiltonian Hy in relative coordinates. In- 
troduce {vi, . . . ,vn) g given by 

1 ^ 

Vi -.^ q^+i - {1 < i < N - 1) and ■= —^qi. 
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3 RELATIVE COORDINATES 



Then v = Mq is a linear change of the coordinates qi, 
matrix 

/ -1 1 ... \ 



, qpf with the N x N- 



M 





.. 





-1 1 

... iv-iy 



The variables u — (ui, . . . ,un) conjugate to v — {vi, 
hy u = {M^)~^p. {M'^)~^ can be computed to be 



, Wat) are then given 



(MJ')-i = _ 



/ 1 

2 



V7V ., 



/I 



1 1 



Nj 



\o... 







1 
...0/ 



(8) 



We already have mentioned that the total momentum J2j=i Pj = is con- 
served, and we write its constant value as N-P. The motion of the center of mass 
Qj — ^Af is linear and therefore unbounded. It turns out that Hy can 
be expressed as a function of the canonical variables {v,u) = {vk,Uk)i<k<N-i 
and P. To express Hy in terms of {v,u), note that by ([5]), Uk = kP — Y^j^iVj 
for any 1 < fc < — 1. Hence 

Pi = -Ml + P; Ptv = Mw-i +P;Pk^ {uk-i -Uk) + P (2 < fc < iV - 1) 
and thus 



N 



NP^ 



1 



= —j^ 1- 2 ("l + ^ "2)^ + . . . + (mW-2 - MAT-l)^ + u^_i) 



Moreover, using that q^+i - qp^ = q^ - q^ = - Ylk^ii'lk+i - Ik) one gets for 
any s G Z>i 



N 



N-l / N-1 N 



fc=i 



fc=i 



Af-l 

fc=i 



(JV-l N 
J2vk 
k=i J 



Combining the two expressions displayed above yields 

Hy — — ;:; V Hy , 



where 



N-2 



Hv^\[ul+Y.^ui+i-uif+ul_A + M E ^fc 



1=1 



k=l 



E 

\k=l / 



(9) 
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N-1 /N-1 



k=l \k=l 



Note that for any values of P, a, and /3, the point {v,u) = (0,0) is a critical 
point of the Hamiltonian Hy . We will see in the next section that it is an elliptic 
fixed point. 

4 Linearized BirkhofF coordinates 

We now compute the Birkhoff normal form of the FPU Hamiltonian Hy up to 
order 4 near the fixed point {v, u) = (0, 0), taking the expansion ([9]) of Hy as a 
starting point. Write Hy as 

Hy = Hu+H,, (10) 

where and Hy denote the u- and w-dependcnt parts of ^ , respectively. Note 
that the Taylor expansion of Hy at (v,u) = (0,0) is not in Birkhoff normal 
form up to order 2. In a first step we therefore want to choose a linear canonical 
transformation {^k,Vk)i<k<N-i '-^ {vk,Uk)i<k<N-i so that when expressed in 
the new variables (^,?7) = {£,k,'nk)i<k<N-i, the FPU Hamiltonian is in Birkhoff 
normal form up to order 2. 

The proposed canonical transformation has naturally come up in our earlier 
work on the Toda lattice f5^. It turns out that (almost) the same transformation 
works for any FPU chain. 

It is convenient to use complex notation for S,k, rjk {I < k < N — 1), 

Cfe ~ '^-'^ "Tf ^ 

The minus sign in the definition of (t is chosen so that d(k A dC-fc = id£,k A drjk- 
The vector ^ = (Cfc)i<|fc[<Af-i is an element in the space 

Z:={z=izk)i<ik\<N-ieC^^-^ -.z-k^z^ yi<k<N-l}. (11) 
The components of C satisfy the identity 

e-.'^/A^C. + e-^-^'/^C-k = V2 (cos Ck + sin (^^^ . (12) 

For the rest of this paper, we use the notation 

Afci-lsin^^l^ (0< |fc| <7V-1). (13) 

The proposed transformation Z — > M^^^'^, C 1-^ (u, u) is given by the formu- 



las 



= ^ E ^'^Cfc, (14) 

V-i' t ^ \ \ ^ AT 1 



l<\k\<N-l 
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4 LINEARIZED BIRKHOFF COORDINATES 



ui+,iC)-n{C) = ^ Afee2-*"=/% {l<l<N-2), (15) 

^ l<\k\<N-l 

-u^^,iC) = ^ Y: A.e-^(^-i)'=/^a. (16) 



l<|fc|<A'-l 

and 

= ^ E Afce2""=/^e-*-'=/^a (l</<iV-l). (17) 

^ l<|fe|<JV-l 

Note that p6|) is actually a consequence of p4)l and ([15]): The left and right 
hand sides of p^ - p^ both add up to 0. From the fact that the Toda lattice 
is integrable it follows that the transformation above is a canonical linear iso- 
morphism [?]. To make this paper self-contained we directly verify that this is 
indeed the case. 

Lemma 4.1. The linear transformation Z M^^^-^, C {v,u), as defined by 
|J4[)-(|77|), is bijective and canonical. 

Proof. First let us show 

{viiC),U^{C)}^tSlm, (18) 

{viiC),vUC)^0, (19) 
{ui{C),Um{C)} = (20) 

for any 1 < l,m < N — 1. Since {v,u) are canonical coordinates on M-^(^~^', 
the proof of (|18p amounts to showing that 

^ f dvi du„i dvi dum \ _ ■ ^ 
for any 1 < m < - 1. It follows from (I14l)-((T7l) that for any \ <k < N 
dCk VN ' dC-k ^ 



o , m — 1 p. , m — 1 

^'^rn _ M yr-^ 27rijk/N ^^ra _ ST^ -2'!zijk/N 



Hence 

dvi dum dvi dur, 



dCk dC^k dC-k dCk 



• 2 / m — 1 m — 1 

Afc I g7ri(2/-l)fe/Af ^ g-27riifc/JV _ g-7ri(2/-l)fe/JV ^ ^2Trijk/N 



N 



3=0 j=0 



i f 2fc7r(l- (/- 7)) 2kn(l-j) 

^ ' cos ^ — - COS ^ ' 



E 



N ^ \ N N 

j=o 

where for the latter identity we used that 2 sin x sin y = cos(a; — y) — cos(a; + y) . 
Taking the sum over k and changing the order of summation then leads to 

E^'^ f dvi dura dvi du„i\ i '^^^Y 2fc7r(l-(?-j)) 2fc7r(/-j) 

. m — 1 
I 



N 



J2 ^('^'-jm - ^1-1,0) 



m~l 

= i ^ - Sij) = i{dim - Sio) = iSim, 

j=o 

as claimed. To prove and (PO)) one argues in a similar way. From 

||^D|) it immediately follows that the linear map ^ 1-^ {v, u) is bijective and 

canonical. □ 

We now compute Hy in terms of the new variables C,. According to the 
decomposition (fTO|) . we compute Hu{C) and i?t;(C) separately. To obtain Hu{C), 
we substitute into the expression i (uf + ~ + ""at-i 

and get 

^"(o^^eY E A,e--/-a- 



2N 

1=0 \l<|fc|<A'-l 



m-1 



i<\k\^k'\<N-l \l=0 / 



k'- 



Using again that Eilo^ e^'^'"^/^ = Ndko and Afe = A_fe for any < < iV- 1, 
one obtains 

Huic) = E AfcCfeC-fc- 

k=l 

Before computing Hy(Q, we simplify its expansion in terms of the variables 
{vk)i<k<N-i- Define vq by the expression on the right hand side of (jl7p evalu- 



10 



4 LINEARIZED BIRKHOFF COORDINATES 



ated at / = 0. Note that 

Af-l , /N-l 



1=0 " l<|fc|<Af-l \ (=0 



Hence X^zl^i^ — and therefore 



(=0 



Substituting the expression p7p for in the quadratic term in the expansion 
([2T|). we get 



Af-l /Af-1 



2 ^ ' 2N 

1=0 l<|fc|,|fe'|<Ar-l \l=0 

N-l 



A:=l 

where we again used that = X^k and ^^^^ gZTriifc/iv _ A^(5;.q for any < 
\k\ <N-1. 

The terms of third and fourth order in Hy are treated similarly. Combining 
the above computations leads to 

HviO = G2 + aGs + /3G4 + O(C') 

with 

N-l 

G2:=2E A^aC^fe, (22) 
fe=i 

G3:=--i= 5] (-l)('+'='+'=")/^AfeAfc,Afc.aa'a", 

{k,k',k")£K3 

^^■■^^ E (-i)''+''^"^'"''/'^A,Afe,Afe.Afe.<fea.a.a,., (23) 



247V 

(k,k\k",k"')£K4 



where 



Ks {(fc, k\ k") cZ^ : 1 < |fc'|, \k"\ <N ~l (24) 

and k + k' + k" = Q mod iV} 

and 

Ki:^{{k,k',k",k"') eZ*: 1 < |/c|,|fc'|,|fc"|,|fc"'| < TV - 1 (25) 

and k + k' + k" + k'" = mod N}. 

Note that G2, G3, and G4 are independent of a and So indeed, Hy is in 
Birkhoff normal form up to order 2, and it follows that C = is an elliptic fixed 
point of the corresponding Hamiltonian system. 
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5 Proof of Theorem 11.11 

We now begin by transforming the Hv{C) ij^to its BirkhofF normal form up to 
order 4. Here we follow a standard procedure - see e.g. section 14 in [7]. The 
phase space Z, defined in (fTTjl . is endowed with the Poisson bracket 

l<|fc|<JV-l SK S K 

where ak = sgn (k) is the sign of k. The Hamiltonian vector field Xp associated 
to the Hamiltonian F is then given by Xp — iJ2i<\k\<N-i o"fc ^^^^ ■ With a 
first canonical transformation we want to eliminate the third order term aG^ in 
HviC)- By a by now standard precedure we construct such a canonical trans- 
formation on the phase space Z as the time-l-map ^I'l :— X^p \t=i of the flow 
■^aFs of a real analytic Hamiltonian aF^ which is a homogeneous polynomial in 
Cfc ^ \k\ < N — 1) of degree 3 and solves the homological equation 

{G2,aF3} + aG3 = 0. (26) 

To simplify notation we momentarily write F instead of aF^ and H instead of 
Hv- Assuming for the moment that ((26)) can be solved and that Xp is defined 
for < t < 1 in some neighbourhood of the origin in Z, we can use Taylor's 
formula to expand H o Xp around t = 0, 

r* d 

HoXi = HoXl+ / —{HoX%)ds 
Jo ds 

= H+ I {H,F}oX^pds 
Jo 

= H + t{H,F}+ ds ds' — {{H,F}oXi) 
Jo Jo ds' 

+ t{H,F}+ [ {t- s){{H,F},F}oX'pds. (27) 
Jo 

When evaluating this expression at i = 1, one gets 

iJo5'i==G2 + {G2,F}+ / {l-t){{G2,F},F}oX'pdt 

Jo 

+aG3+ [ {aG3:F}oX*pdt + l3Gi + OiC^). 
Jo 

Using that {G2, F} + G3 = 0, the latter expression simplifies and we get 

i7 o *i = G2 + / t {aGs, F} oX*pdt + (3G4 + O(C^). 
Jo 

Integrating by parts once more and taking into account that F = aF^ is homo- 
geneous of degree 3 one obtains, in view of (|77)l. 

iJy o *i = G2 + i{aG3, aFa} + PG^ + O(C'). (28) 
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5 PROOF OF THEOREM ?? 



Note that {G3, ^3} is homogeneous of order 4. Hence our first step is achieved. 
It remains to solve ((26)) . Since G3 contains only monomials with {k,k', k") G K3 
(cf also i^s need only contain such monomials, 

^.(3) 



which leads to 



^3 - ^ Pkk'k"CkCk'Ck" 
{k,k',k")GK3 



dFs 



where 



l<|fc|<7V-l ^ 

= ^ (.sfc + sfc, + sfc.)^ifc'fc"aa'a", (29) 

{k,k',k")£K3 



Sk 2a-feAfc = 2sin^. 



The following result is due to Beukers and Rink (cf. [HUE]): 
Lemma 5.1. (IJll \J^) For any (fc, k', k") € K^, 

Sk + Sk' + Sk" ^ 0. 



Let us remark that Lemma lOI also follows from the integr ability of the Toda 
lattice (cf. [5]). We include the self-contained proof due to Beukers and Rink. 

Proof. Suppose that (k, k' , k") £ Kj, satisfies Sk + Sk' + Sfe" = 0. It follows from 
fc + fc' + fc" = mod N that either Sk" — —Sk+k' or Sk" = Sk+k' , according to 
whether k + k' + k" = or k + k' + k" = N mod 27V. 
In the first case, it follows that 

2z sm — + 2i sm — - 2t sm + — j = 0. (30) 

Setting X := and y := e^v^, one can rewrite ([30|) as 

= x - - + 2/ - - - + — = (1 - .x)(l - ?/)(l - xy) — . (31) 
X y xy xy 

It follows that any solution of (pij) contradicts the assumption 1 < |fc|, |fc"| < 
A^ — 1. Indeed, solutions with x = 1 (i.e. fc = mod 2N), y — I (i.e. k' = mod 
2iV), or = 1 (i.e. k + k' = mod 2iV and thus k" = mod 2A^), contradict 
this assumption. 

In the second case, we have instead of (l30l) 



13 



With x,y as above, it now follows from ([5^ that 

= X - - + y - - + xy - — = -{1 + x){l + y){l ~ xy) — . 
X y xy xy 

Again we conclude that any solution of ([5^ contradicts the assumption 1 < 
|^Mfc'|,|fc"| < N - 1. Indeed, solutions with x = -1 (i.e. k = N mod 2N), 
y ^ -I (i.e. k' = N mod 2iV), or xy = 1 (i.e. k + k' = mod 2N and thus 
k" = N mod 2N), contradict this assumption. □ 

By Lemma one can define F3 as follows 

ik,k\k")eK3 

otherwise 



''■^kk'k" 




Then {G2, aF^} + aG^ = 0. Written more explicitly, the nonzero coefficients of 
F3 are 

(3) _ (_i)(fc+fe'+fe")/JV ^Isin^sin-^sin-^l 

" 67^ 2sin^ + 2sin^ + 2sin^' 

In a second step we normalize the 4th order term /3G4 + ^{6*3, F3} in (pS)) . 
We decompose this sum into its contibution to the Birkhoff normal form and 
the rest, to be transformed away in a moment. Let us first compute ^{Gs, F3} 
in a more explicit form: 



a2 p. i 2 dGs dF3 1 2 dG^ 8^ 

2 2 , ^ aCfc oc-fc 2 , oCk oC-k 

i< fc <JV-i K i< fe <Ar-i ^ 



27V ^ 

l<|fc|<Ar-l 



/ 

3 




3 , AfcA/'A„ 



r+r' 5fcA;AmA/' A„i' 



l<|fe|<W-l l<|!|.|m|.|l'l,l™'l<iV-l 
I' +m' -r' N = k 

where for the latter equality we used that 2afcA^ = Sk- Setting 
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one then gets 



„2 



l<|fcl<Ar-l ' + = ' W I m ;/ K 

~ + = mod iV 

2 JV-1 \ \ \ \ 



16N ^ -1 + (sj' + Sm')/sfe 

l^ + 7n' = k mod N 

k = l l + m = k mod N \ ' ' III // K. 

I' + m' = -~k mod N 

"EE'' 



felt i^mJ^modN^^^ + (sr+SmO/sfc -1 - (S;+Sm)/Sfc 

l'+m' = k mod N 

•(— l)^''"''"'A/A„lA^'A„^'C^Cn^C^'Cm'• 
Note that for k = I' + in' + r'N with 1 < A; < TV - 1 and r' e Z we have 

Sfc = |s;'+m'|- 

Introducfjl for any (Z, m, Z', to') G 

1 1 



, if Z + TO ^ modiV 

Clml'm' — \ (.04; 



otherwise. 



We then get 



2 

{G3,i^3}= ^ Q„z-™-(-l)^""''"'AiA„ArA™<,CmO'Cm'- (35) 



2 ' ' ^ 16iV 

(i,m,i'.m')eJf4 



Combined with formula (|23l) for 6*4, the quantity /3G4 + ^{Ga, i^a} becomes 
1^ ^ (-l)"""""''"'(/3+^Cfefc-fe-a"')AfcAfc.Afc..Afc..<fea'a"Cfe"'- (36) 



24iV ^ ^ ' ' 2 



We now decompose ([36)) into its contribution to the Birkhoff normal form of Hy 
and the rest, and we denote by ttn the projection onto the former one, whereas 
the latter one will be (partially) transformed away by a second transformation 

^'2. 



■^To keep the formula for Ci^i/^i as simple as possible we have not symmetrized the coef- 
ficients Ci^lim'- 
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2 

„2 



Lemma 5.2. The normal form part of (3Gi + ^{G^, Fy,} is given by 

2 

7rAr(/3G4 + ^{G3,F3}^ 



Proof. The indices A:, fc', k" , k'" of the terms in f3Gi + ^{Ga, F3} contributing 
to the normal form satisfy {k,k' , fc", k'") E , where 

{(fc, k', k", k'") eKi\3l<l <m<N -I such that 

{k, k', k", k'"} = {I, ~l, m, -m}}. (38) 

In the case I — m, {I, —I, I, —1} in p8)) is viewed as a set-Uke object whose two 
elements I and — / each have multiplicity two. 

We investigate n^f^Gi) and ttn{^{G3, F:}}) separately. Let us start with G4. 
We distinguish the cases I — m and I ^ m in . For I — to, there are (2) — 6 
distinct permutations of {k, k' , k" , k'") in , whereas for I ^ to, all 4! = 24 
permutations of m, — Z, — m) are distinct. Hence we have 

n i 

^^(/3G4) = -^ 6 5]Af|0|' + 24 AfA^lOnCn 



;=1 l<i<m<A'-l 
' N-1 



= ^lE^'lCd' + 2 AfA^IOnUl^l- (39) 

y i=l l<i5^m<A'-l 

Now let us compute '!^N{\{Gz,F-iy). We have to single out the matches of 
(|38p for which in addition the coefficient Ckk'k"k"' in (|35p does not vanish, i.e. 

fc + fc' ^ mod iV and k + k' + k" + fc'" ee mod N. 

Hence there are two quadruples (fc, k' , k" , k'") in K^^ which satisfy these addi- 
tional conditions, 

+ fc" = k + k'" = , > 

k' + k'" = °^ fc' + fc" = ■ ^ ' 

In both cases, we have st" + Sk'" = —{sk + Sk' ), and therefore p4|l reduces to 

-2\sk+k'\ 

Ckk'k"k"' = -. r— • (41) 

\Sk+k' \ + Sk + Sk' 

Note that ((4T|) remains valid for k + k' = N, since in this case Sk^k' — and 
Sk + Sk' > as fc and k' must satisfy l<fc,fc'<iV— 1, but not for k + k' = 0, 
since in this case \sk+k' \ + Sk + Sk' = 0. 
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We first compute the diagonal part of ttn{^{G3, -F3}). In this case, the two 
possibihties in (|40|) coincide and the solutions are 

(^'^'^"'^"')-{H,-i,l,Sf ' (42) 

where 1 < ^ < — 1. The sum of the coefficients Ckk'k"k"' for the two cases 
hsted in (US) is 



\\s2i\ + 2si \s2i\-2sij s^i-Asf N 

We now turn to the off-diagonal part of ttn[^{G3, F^}) . The quadruples 
(fc, k', k", k'") e Ki satisfying (gO]) for given {I, m} C {1, . . . , TV- 1} with I < m, 
{k,k') = (±;,±m), and {k" ,k"') = (±;,±to), are 

{{I, TO, — —to) 
(— Z, —to, ^, to) 

The remaining matches are obtained from (I43p by permuting the first and second 
or the third and fourth columns on the right hand side of (1431) . bringing the 
total number of all matches to 16 = 4 • 4. Note that by formula (|4T|) . these 
permutations leave the value of the coefficients Ckk'k"k"' invariant. Taking the 
sum of the coefficients Ckk'k"k"' for all the quadruples listed in we obtain 

'^{ci^rn, — l, — m ^ ^l, — m. — l,m ^" ^—l.7n,l,—7n 4" ^—l,—m,l,m) 
„( \si+m\ , |si-m| 



\si+m\ + Si + S„i \si-„i\ + Si - Sr. 
\si-m\ , \si+m\ 



\si-m\~Si+Sm |s;+m|-S;-S 

= -16 



l-m - (si - S™)2 sf^^ - (S/ + SrnY 
_ — 16(23^_^„3^^„^ — S^_„^(5; + Sm)'^ — sf^„^{si — 5m)^) 

sLi^H™ + {Sl-Sm)^{si + Sm)^-Sf_^{si + Sm)^-Sl^{si-Sm)'^ 

= -16, 

since s^_„jS^_^„j = (s/ — s„i)^(s/ + Sm)^. Collecting terms, we thus have 



-a^(y{G3,F3})=^ I ^4cos2 1^101^-16 E AfA^IOnC. 

^ ' 1=1 l<l<m<N-l 



nl. 



a" 



N-l 

^1 E(i-^')ioi'-2 E A?A^ionc™p|.(44) 

; = 1 l<i7tm<Af~l 
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Adding up ^ and dH]), we obtain □ 

Now we want to remove [as much as possible of] the term (Id — 7rjv)(/?G4 + 
^{Ga, F^}) from the Hamihonian Hyo'i'i given by by a second coordinate 
transformation ^'2. In view of formulas and (pS)) for G4 and ^{Ga, F3}, re- 
spectively, and in complete analogy to the first step we look for a transformation 
^2 of the form ^2 = Xpjt=i with 

Fi= ^ -F'it'fc"/c"'CfcCfc'Cfc"Cfc"': 

{k,k',k",k"')eK4\K^ 

where F^'^d^k' ^k" ,k"') = ^(k.k' .k" ,k"') permutation a{k,k' ,k" ,k"') of the 

quadruple (fc, /c', A:", /c'") e K4\K^ . We would like to determine the coefficients 
of F4 in such a way that 

{G2, Fi} = -(Id - nN)iPGi + ^{Ga, F3}). (45) 
As in (pg)) one gets 

{G2,F4} = -z ^ (sfe+,5fc.+Sfc..+,Sfe..0^ifc'fe"fe"'aCfe'Cfe"a"', (46) 

(fe,fc',fc",fc"')Gi<'4\if|' 

and equation ([l5|) combined with ([55]) leads to 

,(sfe + s,,+sfe„+sfe.0^^ifc^fe"fc'" (47) 

for any quadruple (fc, A:', /c", fc'") in \ . Here Ckk'k"k"' denotes the sym- 
metrized version of Ckk'k"k"', 

Ckk'k"k"' ^ X] Ccr(*:,fc',fc",fc"')- (48) 

The following lemma due to Beukers and Rink (cf |llj ) determines the quadru- 
ples {k,k' ,k" ,k"') e if 4 \ for which Sk + s^/ -I- Sk" + Sk"' = 0. Let us 
introduce 

Kr K+, U C 

where 

f TV 
if^t. — {{k, k', k", k'") e K4\3l eN -.1 <l < — so that 

N N 1 

{k, k\ k", k'"} = {±1, ±l^N,-Tl,-Y^l}}- 

Note that if N is odd, then Kres — 0- 
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Lemma 5.3. (^]) Let (fci, /fc2, ^3, ^4) € \ . Then 

Sk + Sfc' + Sk" + Sk'" = if and only if (k, k' , k" , k'") £ K^'^ 
In particular, if N is odd, then + s^,/ + sj.// + Sk"> ^ 0. 



For the convenience of the reader a detailed proof of Lemma 15.31 is given in 
Appendix [A) 

By Lemma 15.31 if N is odd, (|T7)) can be solved for any 
K4 \ determining the coefficients i^(fc,fc',fc",fc"') with {k,k' , k" , k'") E K4 \ 

in such a way that Fj^^i i^t k'") ~ ^{k\' k" k"') ^'^^ permutation 
C7(A;, k', k", k'") of (fc, k' , k", k'") e K4 \ K^. With this choice of F4 the canon- 
ical transformation ^2 is then defined by Xpjt=i- Composing and ^2, we 
obtain the transformation S := o vj/j. We have proved the following 

Proposition 5.4. Assume that N > S is odd. The real analytic symplectic 
coordinate transformation C, — defined in a neighborhood of the origin in 

Z, transforms the Hamiltonian Hy into its Birkhoff normal form up to order 4. 
More precisely, 

Hvo^ = G2+^N (^/3G4 + y {G3, Fa}^ + ©(z^), (49) 

with G2 and TiNiPG^ + ^{6*3,^3}) given by 12^) and IS7\ I, respectively. 
Theorem 11.11 can now be proved easily. 



Proof of Theorem Proposition 15.41 provides the Taylor series expansion of 
Hv in terms of the actions 

I ^ {Ik)i<k<N-u 4 = ^^^. (50) 
More precisely, fly oE = Ha,i3{I) + 0{z^), where H^.^il) is defined by 



k=l k=l i^^m 

l<l,m<N -1 

and Afc = | sin 7^ | ^ . This proves Theorem 11.11 □ 



6 Proof of Theorem 11.41 

Now we assume that N is even. To obtain the normal form of the FPU 
Hamiltonian as claimed in Theorem 11.41 we continue the investigations of the 
previous section. According to Lemma 15.31 equation (|17|) might have no so- 
lution F^f^if^uj,,,, for {k,k',k",k'") e K^'^^. We first compute the projection 
■KresifiGi -\- ^{Gsj^s}) of (HG^ + ^{03,^3} outo those terms which are in- 
dexed by quadruples (fc, k' , k" , k'") G K^^'^, i.e. the projection onto the resonant 
non-normal form part of /3G4 + ^{63, F3}. 
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Lemma 6.1. Assume that N is even. The resonant non-normal form part of 
/3G4 + ^{GsjFa} is given by 

TTres (pG, + ^{Gs,F,}] = -^^{R + R.) (52) 



where 



R := 

i<;<4 



with Cfc :— 2 cos ^ for 1 < k < N — 1, and 



R.^{ ^ ' ' ' (54) 
otherwise. 

Proof. Consider the formula ((36|) for /3G4 + ^{63,^3}. At this point we need 
to consider the symmetrized version (|48l) of the coefficients Ckik'v defined by 
We claim that for any (fci, A;2, ^3, ^4) S K'^'^'' 

^ 2 

^ ^ 5Z '^*-V(l)fe.T(2)fe,T(3)fe,T(4) — (55) 
■ creS4 



^kik2k3ki 



Observe that Ckik2k3ki is invariant under the transpositions ki k2 and ^3 
^4. Hence ([55]) follows once we prove that 



4 {Ckik2k3k4+Ckik3k2k4 +Cfcifc4fc2fc3 +Cfe2fc4feife3 +Cfe2fe3felfc4 + Cfcg ^4 fej fea ) — " 16 . (56) 

Note that any element (fci, k2, k^, k4) £ K^"^^ is, mod 2N , a permutation of an 
element of the form {I, -N + Z, N/2 + Z, -iV/2 + Z) with 1 < |/| < N/A. For such 
quadruples one gets by a straightforward computation 

Ckik2k3k4 + Ck3kikik2 = —2 — 2 — —A 

and, with c; = 2 cos ^ , 

^kiksk2k4 ~l~ Ck2k4kik3 



2+{si+ci) 2-{.si+ci) S2i 
as well as 

4 4 8 

Ckik4k2k3 ~^ Ck2kskik4 ^ ~^ 



2 + (s;-q) 2-(s;-q) S2i' 



Substituting these three identities into the left hand side of ((56l) leads to the 
claimed identity ([56]) . 
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Moreover, by the definition ([55]) of Simi'm' one lias for any (fci, fc2, ^3, k^) G 
Kl"'' and any a S 5*4, that £k„^r^k^^2)k^^3)k^^i~, = ±1 and hence 



Further, 



-1 


^'=,7(1) '=,7(2) '=(7(3) 




= -1. 




. Itt In 


1 


21tt 














~ 2 





Combining all these computations we get 

„2 



(ki,k2,k3,ki)£Kl''^ 



247V 



247V 



247V 



\ 4 4 4 4 / 



(57) 



only if -Jgn 



4A^ 



{R + Rn), 



with i? and i?jv as defined by ((53l) and ((54l) . respectively. Hence Lemma [6. II is 

4 

proved. □ 



By Lemma 16.11 if N is even, equation (j47p can be solved for any quadru- 
ple {k,k',k",k"') (z K4,\ {K^ U Kl''') in such a way that Flf^ ,^, ^^„^^,„^ = 

P{kk' k" k'") foi' any permutation a{k, k' , k" , k'") of {k, k' , k" , k'"). With this 
choice of F4 the canonical transformation ^'2 is then defined by Xp \t=i- Com- 
posing ^'i and 5*2, we obtain the transformation S := o v[/2 and have proved 
the following 

Proposition 6.2. Assume that N is even. The real analytic symplectic coor- 
dinate transformation ^ = S(z), defined locally in a neighborhood of the origin 
z = in Z, transforms the Hamiltonian Hy into the resonant Birkhoff normal 
form up to order 4, 



a 



i/y o S = G2 + TTJV /3G4 — {Ga, F^} + nres PG^ + — {G3, F3} + O(0^), 



with G2, iTN{l3Gi + ^{G3,i^3}), and -KresiPG^ + ^{G3,i^3}) given by 
PTp , and h52\l . respectively. 
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Proof of Theorem \1.4\ We start with the formula for Hy o S given by Proposi- 
tion [621 and treat the normal form terms G2 + irNiPGi + \ {Gz, ^3}) and the 
resonant non-normal form terms TTresiPGii + ^{G^, F3}) separately. With the 
action variables / = (/fc)i<fe<jv-i defined by we see that G2 + ttn{(3G4 + 
^{G3,F3}) = Ha.j3{I), where Ha. j3{I) is defined by ([^ . Concerning 7rres(3G4-|- 
^{Gs, F3}), we first express it in terms of the real variables {xk,yk)i<k<N-i, 
related to the Cfc's by Xk = (Cfc + C-fc)/2 and yk = (C-fc - Ck)/2i. Note that 

CiC-N+i(f_iC_N_i + C-iCN-iC-f+iCf+i 

= 2Re{CiC-N+iCf-iC-»-i) 

= \{ {xiXN-i + ViVN-i) {xN_iXN^i + 2/« _;y« 

- (xiVN-i - XN-m) {xN_iyN^i - XN^iyN_^ ) 
= 2 [JiJn_i - MiMn_^ , (58) 
where for any 1 < fc < — 1 

Jfc := ^ {xkXN-k + ykVN-k) and Mk := ^ [xkVN^k - XN-kVk) ■ 



Hence i?, given by ([53|) . can be expressed in terms of Jk and Mk as follows 
R{J,M) = 52/ (0C-w+;Cf-;C-f-( +C-iCiV-iCf +;C-f +/) 



i<;<f 

2l7r 



= 4 E sin^ (j,J«_, - M,M«_;) . (59) 
i<i<f 

Similarly, if f G N, one concludes from (|58p that i?™., given by ([M)) . can be 
expressed as 

Rn{J,M) = i fciC!3iv +C|vC!«) = Jl-Ml. (60) 



2 

Theorem [Til now follows from the formulas ([SI]), ([Ml), and JMl). □ 



7 Proof of Theorem 11.31 

To analyze the Hessian Qa,p of ([5T|) at / = we repeatedly encounter matrices 
of the form E + diag(/Lti, . . . , fj,N-i), where E is the (TV — 1) x {N — l)-matrix 



E := 



/I 



(61) 



and (/ife)i<fc<jv-i are given complex numbers. The determinant of the matrix 
E + diag(/ii, . . . , fiN-i) can be explicitly computed. 
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Lemma 7.1. Let {lJ-k)i<k<N-i be given nonzero complex numbers. Then 



(N-1 ^ \ Af-1 
1 + X! — 'IT 
k=i '^V fc=i 



Mfc- (62) 



In particular, E + diag{iii, . . . , /ijv-i) is regular if and only i/X^I^i^ 77^ 7^ 

Proof. Expanding dct(_E+diag(/ii, . . . , ^n-i)) with respect to its rows it follows 
that 

N-l N-1 

det(£; + diag(^i,...,/iA,_i)) = H ^'■^ + X! Il'^'- 

k=l k=l l^k 

This leads to formula JMl)- □ 



First let us treat the case a — 0, [3 0. Using the notation introduced in 
section SI one obtains the following proposition. It improves earlier results of 
Rink [11]. 

Proposition 7.2. Let N be odd and assume that a = in Then the 

following holds: 

(i) The Birkhoff normal form of Hy up to order A is given by — \- Hq,p{I) 
where 

AT-l /n-i 
HoAl) = 2 ^ ^'^'^ + ^ E ^tl'k + 2 E I • (63) 

fe=l \k=l l<l^m<N-l 

(a) For any /3 7^ 0, Ho^/3{I) is nondegenerate at I ~ 0. 

Proof. The Birkhoff normal form ((63|) of Hy is given by the formula (|5ip eval- 
uated at a = 0. To investigate the Hessian (5o,/3 of Hq^p{I) at / = we write 

Qo,/3 = ^APA, (64) 

where 

A = diag I sin — ^ ) (65) 

and 

P^2 - (e - Imn-i 



2 



In view of (l62l) it follows that 



^ ' k=\ 



23 



where by Lemma |7. 11 

detP = 2^-1 (1-2(7V-1)) (-1/2)^-1 = (-l)^(2iV-3) 7^0. 

Hence, if /? ^ 0, detQo,/3 7^ 0, and the nondegeneracy of Hq^p^I) at / = 
follows. □ 

Lemma 7.3. If P < 0, then Qo.p has one negative eigenvalue, whereas if (3 > 0, 
then Qo,i3 has N — 2 negative eigenvalues. In particular, for any /3 G R \ {0}, 
Qo,/? 'is indefinite (and i?o,/3 "i-s therefore not convex). 



Proof. We want to use the decomposition of Qo,0 to show that (5o,/3 can 
be deformed continuously to -^P- Consider for < i < 1 

QoAt) - ^ ^ + (1 - w) P (i A + (1 - 1) Id). 

As i A + (1 — Id is positive definite for any < t < 1 and P is regular, Qo.fsit) 
is a symmetric regular (A^ — 1) x (A^ — l)-matrix. For t — 0, Qo.p{0) — ^P, 
whereas for t = 1, Qa.fj{l) — Qo,/3- Therefore, index((5o,/3) (i-e. the number of 
negative eigenvalues of Qo./s) coincides with index(^P). The eigenvalues of P 
are ni = 2N — 3 with multiplicity one and /i2 = ^ 1 with multiplicity N — 2. □ 

We now turn to the case a ^ 0. 

Proposition 7.4. A ssume that o; ^ in 1^) . Then, for oi fixed, det (5a./3 is 
a polynomial in (3 of degree N ~ 1 and has N ~ 1 real zeroes (counted with 
multiplicities) which we list in increasing order and denote by (3k — (^kict) A < 
k < N — 1). They satisfy < (3i < < (^2 < ■ ■ . < (3n-i and contain the 
l(A^ — l)/2j distinct numbers 

21. / . 9 kTT\^^\ ,^ . , . A - 1 



Moreover 









fsin^ ^"l 


1 









for (3<(3i 

index {Qa,[3) = ^ for (3i < (3 < P2 



for (3 > (3n- 



1 



Proof. Fix a G M \ {0} and consider the map (3 ^ dei{Qa.p). It follows from 
(|5T|) that det{Qa,i5) is a polynomial in (3 of degree at most A — 1, 



AT-l 



det(Qa,;3) = ^ q](3\ 



where go = det((5a,o) and qn-i = det((5o,i)- By Proposition l7.2[ det(Qo,i) 7^ 0' 
hence the degree of the polynomial de\.{Qa.p) is A^~ 1. We claim that dei{Qa,p) 
has at least A^— 1 real zeroes (counted with multiplicities). For \(3\ large enough, 
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index (Q a, /s) is equal to index((5o,/3)- By Lemma [7.31 index((3o,/3) is N — 2 for 
/3 > and 1 for /3 < 0. Hence there exists R > such that index{Qa^p) = N —2 
for any fi > R and index((5£j,^) = 1 for any [3 < —R. For /? = a^, Qa,a^ is a 
positive muhiple of the identity matrix, hence index (Qq, ,-,2) — 0. It then foUows 
that mdex{Qa^0) must change at least once in the open interval (—00,0^) and 
at least N — 2 times (counted with multiplicities) in (a^,oo). Since a change 
of index(Qct^^) induces a zero of det{Qa,f3) (counted with multiplicities), our 
consideration shows that /3 1-^ det{Qa,i3) has at least — 1 real zeroes. Thus 
f3 det((5a,/3) has precisely A^ — 1 real zeroes and we have /3i(a) < < f32{ct). 

Next we prove that /3i(a) > 0, i.e. that Qa.p is regular for any (3 < 0. Write 
Qa /3 as a product, 

Qa,l3 = A Pa^fi A, 

where A is given by ([65|l and Paj is given by 



l<k<N-lj 




where E is given by (pT|) and 

7(a,/3):= 



v2 



As -00 < ^ < it follows that < 7(0, ^) < 1 and fik = -5 (l + < 

for any 1 < fc < — 1. Lemma [7.11 savs that Pa^p is regular if f{j{a,(3)) ^ 
where 



KTT \ 



/(7) :-l-2E + 



in the interval < 7 < 1. Note that 7(7) is increasing in < 7 < 1 and /(I) 
can be estimated as follows. Using that N is assumed to be odd one has 



sm^ ^ sin^ (^'^^'^ 



f fl) = 1 - 4 V — < 1 - 4- 



fe=i — sm ^ 1 + sin^ ^ 

cos' W 



1-4- ^ = -3 



^ 1 + C0S2 1 + C0S2 ^ ■ 



As for > 3 



4 4 5 

1 + COS^ 2^ 1 + COs2 I 7 

we conclude that /(I) < 0. Hence we have shown that 7(7) < for < 7 < 1, 
and therefore Pa,p is regular for /3 < by Lemma [7.11 Hence we have proved 
that < /3i(a). 
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Finally introduce fik '■= ^5(1 + 7(0^, /?)/ sin^ ^) and note that for /? with 
7(0;, /?) = — sin^ for some 1 < fco < one has fiko = fJ-N-kg = 0- As 

ko ^ N — ko a 1 < ko < \_{N — l)/2j it then follows that Pa. (3 has two equal 
rows and is therefore singular. Note that 7(0;, /3) = — sin^ corresponds 
to /3 = a^{l + sin~^ and we have proved that /3 1-^ det{Qa.i3) has at 
least l(A^— 1)/2-i different zeroes in the interval {a^,oo). The statement about 
mdeii{Qa.f3) easily follows from the above analysis. □ 

Proof of Theorem \1.S[ Part (i) is proved by Proposition l7.41 whereas (ii) follows 
from Proposition 17.21 and Lemma l7.3l □ 



A Proof of Lemma 15.31 

For the convenience of the reader, we provide a detailed proof of Lemma 15.31 
in this appendix. This lemma and its proof are due to Beukers and Rink - see 
([TT]. Appendix A). Recall that K4\K4 C denotes the subset of quadruples 
{ki,k2,k3,k4) satisfying 1 < \ki\ < iV - 1 (1 < z < 4) and ki + k2 + ks + k4 = 
mod N so that there are no integers I, m with {I, m, —I, —m} — {ki, k2, k^, k4}, 
and 

Kr KL U K-,, C K4 

where 

K^^^ := |(A:i,A:2,fc3,fc4) e K4\ 31 e N : 1 < I < — so that 

N N 1 

{fcl, k2,k3,k4} = {±1, ±lTN, — Tl,-Y'f ^}|- 

Note that K^'^'' = if iV is odd. Let us restate Lemma [5731 as follows: 

Lemma A.l. ffllf ) Let (fci, A:2, ^3, ^4) be an element of K4 \ K4 . Then 
(fci, k2, ^3, k4) € K\'^^ if and only if 

kit: fc27r fcsTT fc47r 

sm— +sm— +sm— +sm— =0. 

Let us make a few preparations for the proof of Lemma lA.ll By a straight- 
forward computation one sees that the "only if" -part of the claimed equivalence 
holds: 

Lemma A. 2. For any {ki, k2, k^, k4) € K^'^'^ , one has ^fr — 0- 

So it remains to prove the converse. First we consider some special cases. 

Lemma A.3. Let (fci, fe, ^^3, fc4) G ^4 \ (K^ U Kl^"). If there exist I, 
such that 

(i) {ki,k2,k3,k4} = {l,-l,m,n}, or 

(ii) {ki, k2, k^, ^4} — {I, N — I, m, n} with I < I < N — 1, or 
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(in) {ki,k2,h,k4} = {l,-N - l,m,n} with -{N - 1) < I < -1, 
then 

4 , 

sm— ^0. 

Proof. In case (i), it follows that m + n = N (and thus 1 < m, n < — 1) or 
m + n = —N (and thus —{N — 1) < m,n < —1). Hence in both cases, sin ^ 

and sin ^ have the same sign and ^in = sin ^ + sin ^ 7^ 0. In the 

case (ii), by assumption, m + n = mod N. The case m + n = has already 
been treated under (i). U m + n = N, then sin ^ > for any 1 < i < 4. If 
m + n = —N, then m < 0, and m ^ {—I, -N + I}. Thus n = —N — to < 

and therefore X]i=i ^W' ~ ^ ^^^^ W ^ ^^^^ ^^"n''^ 7^ 0- The case (iii) is treated 
similarly as (ii). □ 

Another special case in treated in the following lemma. 

Lemma A. 4. Assume that (fci, fc2, fcs, ^4) G \ satisfies 

k, + kj ^0 mod N V 1 < i, j < 4. (66) 

// there exist l,n g {ki, k2, ks, k^} with 

In rm 

sin ^ + sm ^ = 0, (67) 



then 



^sin^^O (68) 



implies that (fci, fc2, /ca, fc4) G i^l'^'*. 

Proof. From the assumptions ([55)) - (l67p it follows that there exists 1 < Z < 
iV — 1 so that {fci, A:2, /ca, fc4} = {Z, — iV + Z, to, n} for some m,n G Z. Then 
sin^ + sint-(^ilZL = Q and hence by (UHl), sin + sin ^ = 0. W.l.o.g. 
assume that 1 < to < iV — 1. Then either n — —to or n = —N + to. If n = —to, 
then (/ci, /c2, ^3, A:4) G K^^^ by Lemma [A.BI fiV If n = —N + to, then one has 

'h^2l-N + 2m-N = 2{l + m)-2N. 



Note that 2(Z + to) — 2N cannot be an even multiple of iV, as otherwise Z + to = 
mod N , violating ((66)) . If, in addition, N is odd, then 2{l + to) — 2A^ cannot be 
odd multiple of N. Hence in the case N is odd we conclude that X]i=i ^ 
mod N , contradicting the assumption (fci, k2, fcs, k4) G -fi'4. 

If is even, it is however possible that 2{l + to) — 2N equals ±iV: If 
2{l + m) ~2N ^ N, i.e. I + m = ^N, it follows that f < Z, to < iV - 1, and 
(fci, fc2, /C3, ^4) G i^-^s with I' ^l-f. If 2(/ + to) - 2iV = -N, i.e. Z + to = f , 
it follows similarly that (fci, A;2, ^3, ^4) G ^its with /' = I. In both cases, we 
conclude that (/ci, /c2, fcs, fc4) G X^'^'*. □ 
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In view of Lemma IA.3I and Lemma IA.4I in order to prove Lemma lA.ll it 
remains to show the following 

Lemma A. 5. Assume that (fci, fc2, fca, /C4) G K/^ satisfies I166\} . If for any 1 < 

i,j < 4 

sin^+sin^^O. (69) 

(and thus (fci, fe, ^3, ^4) ^ K4 U K^'^^), then 

4 

sin 

N 



1=1 



To prove Lemma [A. 5 1 let us first rewrite using Euler's formula for the 
sine function, 

E = (70) 

i<lil<4 

where C±j = -j-giifciTr/Af 2A^'th roots of unity. Note that for any quadruple 
{ki,k2,k3, € ^4 \ satisfying ((69l) one has 

+ VI < bi < i/i <4 

Indeed for any 1 < |j| < < 4 one has Im Q+ Im (ji = sin -^^^^ + sin ^-^^^ 
which does not vanish by assumption (|69p . 

Let us first discuss equation (|70|) and its solutions in general. For convenience 
let us rewrite ([70]) as 

Ci + . . . + Cs = 0. (71) 

We are interested in the solutions {Ci)i<i<s of the equation ((7T|) on the unit 
circle S'l := {z G C||z| = 1}. 

We need an auxiliary result which we want to discuss first. Let n > 2 
be arbitrary and assume that the sequence (Ci)i<i<ra Q has no vanishing 
subsums (i.e. J2ieJ for any 7^ J C {1, ... , n}) and satisfies the equation 

n 

5^0 = 0. (72) 

1=1 

Let M G N be the smallest integer with the property that (Ci/Cj)'^ = 1 for all 
^ < i,j < n. Then there exists ^ G 5'^ so that = for any 1 < i < n. 
W.l.o.g. we can assume that ^ = 1. Furthermore, let p'^ be a prime power 
dividing M so that M/p'^ and p are relatively prime and define 

M' AI/p and 77 e^^^/p' . (73) 

Then for any I < I < n there exists a unique integer < /i(/) < p — 1 such that 
Ci^Ci- ?7^^'^ where is an element of the field K Q(e2'^'/^^'). (As Ct' = 1 
there exists < ri < M ~ 1 with Q ~ e^^' . If r; = mod p choose /i(/) = 0. 
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If n ^ mod p choose 1 < < p — 1 so that r; = mod p.) Hence (17^ 

can be written as 

n n p— 1 / \ 

o = Eo = EO'?'''^ = E E cAv' (74) 

1=1 1=1 s=n yiep-i(s) J 

We need the following algebraic fact (see e.g. [13], §60-61): 

Proposition A. 6. The minimal polynomial of rj — e^'^^lv over the field K = 
(Q(e2WA^') is given byXP-r]P if k > 2 and XP-^ + X''-^ + . . . + X + 1 if k = 1. 

We now claim that M is square- free, or equivalently that for any prime power 
p'' dividing M, 

k = l. (75) 

Indeed, equation ([74]) shows that the minimal polynomial of C, has degree at 
most p — 1, which by Proposition IA.6I is only satisfied in the case fc = 1. 
Further we claim that there exists cr g C \ {0} so that 

J2 6=0- V0<s<p-1. (76) 

The existence of such a a follows from Proposition IA.6I As fc — 1 by ([75]) , 
the minimal polynomial of ry over K is given by XP~^ + XP^^ + . . . + A + 1. 
Since this is a polynomial of degree p ~ I the polynomial on the right hand 
side of (|74p must be a scalar multiple of the minimal polynomial. Hence all 
the coefficients X]iep-i(s) have the same value cr G C. As J2i£f^-i-(s) ~ '-'"'f 
the additional property cr ^ follows from the assumption that there are no 
vanishing subsums. Hence we can assume w.l.o.g. that cr = 1. 
Next we claim that 

p<n. [11) 

In other words, possible prime factors of M are bounded by the number of 
summands in ([7^ . To prove ([77)l . note that it follows from ([75]) that for any 
< s < p — 1 there exists \ < I < n such that /i(Z) — s, i.e. the map 
fi : {1, ... ,n} ^ {0, ... ,p — 1} is onto. This establishes ([771) . 

The map induces the partition (jj /j,^^(s))o<s<p-i of the positive integer n 
into p summands, 

p-i 

'^ = Ett^"'(^) (^8) 

s=0 

Lemma A. 7. ('ll'l, Appendix A) For any solution {(i,...X8} of (7l\ ) con- 
tained in without vanishing subsums there exists G such that either 

{Ci, ...,Cs} = {-^a, -^a'} U {^7-' 1 1 < J < 6} (79) 



or 

{Ci, . . . , Cs} = {-^a\ -Ca' • f3\-^a' ■P'\l<l<2}U C/3™}, (80) 
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where the quadruple (i,j,k,l) is a permutation o/ (1,2, 3, 4) and 

2-Ki 2-Ki 2Tvi 

a:=e3, p:=e=, "f :— e . 

Proof. By a straightforward computation one verifies that the sets of the form 
([79]) or (jSO]) satisfy (|7ip . It remains to prove that these are the only sohitions 
of ([71]) of this type. 

We classify the solutions of ([71]) according to the possible values of p, which 
we now assume to be the largest prime dividing M. Since n = 8, by (|77p . the 
possible values of p are 2, 3, 5, and 7. li p = 2, then, by (I75p . Af = 2 and 
therefore there exists ^ G S*^ so that = ±C for any 1 < i < n. In this case 
there exists a solution of ([7^ without vanishing subsums only if rt = 2. (In this 
case, they are given by {Ci, C2} — — 1} with ^ £ .) If p = 3, then M — 3 
or M = 3 • 2, and there exists a solution of ([7^ without vanishing subsums 
only if TT, — 3. (In this case, they are given by {Ci:C2,C3} = Ci^i'^i'^^} with 
^ G S^.) li p = 5, then 77 = /3 in (|73p . Up to permutations, there are the 
following three partitions of 8 into 5 summands, (4, 1, 1, 1, 1), (3, 2, 1, 1, 1), and 
(2, 2, 2, 1, 1). In a straightforward way one shows that the partitions (4, 1, 1, 1, 1) 
and (3, 2, 1, 1, 1) and their permutations give rise to solutions of the equation 
([7T|) with vanishing subsums. E.g. the solutions corresponding to (4, 1, 1, 1, 1) 
are given by ^ • (-/3, -/J^, -/J^, /3, /J^, Z?-*) with ^ e 5\ whereas the 
solutions corresponding to (3,2,1,1,1) are ^- (— i, 1, i, — a/3, — /3'^, /?^) 
with ^ G S^. On the other hand the partition (2, 2, 2, 1, 1) leads to the solutions 

(Ci, • • • , Cs) - e(-a, -a^, -a/3, -a^/?, -al3\ -a'p^p^ l3^) 

with ^ G S^. They are the solutions (HO]) with {i,j,k,m) = (1,2,3,4). Permu- 
tations of the partition (2, 2, 2, 1, 1) again lead to solutions of the type ([80]) . but 
with (i, j, fc, to) given by a permutation of (1, 2, 3, 4). 

If p = 7, then 77 = 7 in ([75)1 . Then, up to permutations, (2, 1, 1, 1, 1, 1, 1) is 
the only possible partition of 8 into 7 summands. The partition (2, 1, 1, 1, 1, 1, 1) 
leads to the solutions 

(Ci,---,C8) =£.{-a,-a'^,j,...,j^) 

with ^ G S^, where we used that 1 ~ —a — a^. They are of type ([79)1 . Any 
permutation of (2, 1, 1, 1, 1, 1, 1) leads to the same kind of solutions. □ 

Lemma A. 8. fU 1^. Appendix A) For any solution {(i,...Xs} of con- 
tained in without vanishing subsums of length 2 but having a vanishing sub- 
sum of length 3, 4, or 5, there exist ^, ^' G such that 

{Ci, . . . , Cs} = Ua'lO <1<2}U U'/J^IO < m < 4}, (81) 
where again a = e^'^^l'^ and f3 — e^'^'/^. 

Proof. Again, one verifies by a direct computation that the solutions (ISTj) of 
(|7ip have the desired properties. It remains to prove that they are the only 
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ones. First note that under the hypotheses of the lemma, vanishing subsums 
of length 4 cannot occur, since the latter ones would imply the existence of 
vanishing subsums of length 2, which by assumption is excluded. Hence, in 
order to find solutions of ((72|) for rt = 8 with the desired properties, we have 
to find all solutions of (j72p without vanishing subsums for n = 3 and n — 5. 
Note that by (j77p . p — n for n = 3 or n = 5. By the considerations in the proof 
of Lemma lA.71 the former ones are given by (CijC2jC3) = C(IjCKiCk^) a.nd the 
latter ones by (Ci, . . . , Cs) = /?, /3^ /3^ with ^, ^' e S^. This proves the 
lemma. □ 

We are now ready to prove Lemma lA. 51 

Proof of Lemma \A.5[ We first select from jTS]), dHOl) and ([81]) all the solutions 
(Ci, . . . , Cs) of (|7T|) which are of the form ([55]) (after multiplication by 2i). This 
amounts to selecting the solutions (Ci, . . . , Cs) of ([7T|) having the property that 
{Ci, • . • , Cs} is invariant under the map <—> —C^^ . It requires to choose ^ and 
^' in ([79|) . ([80| . and (|8T]l appropriately. Let us explain this procedure in detail 
for the solutions of type (|79|) . 

First we rewrite the solution (TfO]) . 

(Ci, . . . , Cs) = C ■ (-a, -a^ 7, 7^ 7^ 7'', 7^ 7^) = (e^) 

V /l<fe<8 

where 



(ti,...,t8) = (6,7,12,18,24,30,35,36) and x e M/42Z. (82) 

The required invariance of the set of the Cfc's under the map C, i— » — C^^ is 
equivalent to the invariance of the set of the (tk + x) 's under the map ti~*21 — t 
(mod 42). Since the set ([5^ of the i^'s is invariant under the map 1 1^ —t (mod 
42), {tk + x\l < fc < 8} is invariant under t i~> 21 — < (mod 42), if wc choose 
X :— ^ or = i. Then the equation reads 

IItti 5771 ISttz IQni 237Ti 277ri ni 37ri 

+ e~ + + + + + e~ + e^^ = 0, 

or sin -I + sin + sin + sin ^ = 0. Choosing all arguments in (0, tt), the 
latter identity reads 

TT 37r TT 57r 

sm — h sm sm sm — = 0. (83) 

6 14 14 14 ^ ' 

For the solutions of type ([SO]) , one gets 

. TT . 137r . Ttt . 37r 

sm h sm sm sm — = (84) 

6 30 30 10 ^ ' 



and 

TT TT IItT 
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Let us briefly explain how ([5^ - ([55)) can be obtained. Note that from the 24 
permutations of (1,2,3,4) in ((80|) . there are only six which lead to different 
sets of the Ci's, since interchanging i and j or k and m leaves the set on the 
right hand side of ((80|) invariant. In the resulting six different cases, we again 
write {Ci,...,C8} = ■ {e2'^*4n,...,e2'^*-3l} with U in R/30Z. Then, up to 
translations, there are only two different types of solutions emerging from these 
six cases. With the appropriate choices of ^, one gets the solutions ([M]) and 
(ESI). 

Finally, for the solutions of type ([5T|) . one gets 

sm- -sm- +sm— -sm— = 0. (86) 
z D lU ID 



The procedure to obtain ((86|) is basically the same as in the preceding cases. 
We write ^ as {Ci, . . . , Cs} = ^ • A • ^™|0 < I < 2, < m < 4} and first 
choose A S 5"'^ so that the set {a', A • /3™|0 < I < 2, < m < 4} is symmetric 
with respect to some axis through the origin, and then choose ^ so that this axis 
is the imaginary axis. 

To finish the proof of Lemma lA.SI we show that all the solutions (fci, fc2, fca, ^4) 
of J2i=i sin ^ = obtained in and the additional ones obtained by 

replacing < a; < tt in sinx by tt — a; satisfy ^ mod N and hence are 

not in K4. 

For the solutions obtained in (1551) - ([55)) . ^ is even. Hence if N is odd, then 
there is no quadruple (fci, k2,k^, k4) E K4 such that ((55)) and (IM)) are satisfied. 
This finishes the proof of Lemma [A. 51 in this case. 

For the rest of the proof, we assume that N is even. If iV = 42r for some 
r g N, (HSl) becomes 

7r7r 9r7r (— 3r)7r (— 15r)7r 

sin — — h sin — — h sin — h sin — = 0, 

42r 42r 42r 42r 

and we have 7r + 9r — 3r — 15r = — 2r ^ mod 42r. Hence the corresponding 

quadruple (fci, A;2, fcs, ^4) is not in K^. For the quadruples obtained by replacing 

< a; < TT in sinx by tt — a; in some of the summands in (|83p . the condition 

Si=i ^ mod 42r amounts to 

±7±9±3±15^0 mod 42 (87) 

for any combination of plus and minus signs. The relations (|87|) are easily 
verified. Similarly, one verifies that the quadruples (fci, fc2, /cs, ^4) satisfying 
(|84l) . ([85]) . or ([86]) are not in K4 by showing that 

±5±13±7±9^0, ±5±1±11±3^0, ±15±5±3±9^ mod 30, (88) 

again for any combination of plus and minus signs. Hence we have shown that 
none of the solutions (fci, k2, k^, k4) of ([68)1 is an element of K4. This completes 
the proof of Lemma fA.Sl □ 

Proof of Lemma \A.l\ The claimed statement follows from Lemma |A.2[ IA.3[ 
\EM andEl □ 
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